Discrete Spectrum in Bose Condensation 
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A quantitative description of the whole process of condensation of bosons in an harmonic trap is 
given resorting only to Gibbs and Bose postulates, without assuming equipartition nor continuum 
statistics. Below Tc discrete spectrum theory predicts for the thermo-dynamical variables a behavior 
different from the continuum case. In particular a new critical temperature Td is found where the 
specific heat exhibits a A-like spike. Numerical values of the relevant quantities depend on the 
experimental set-up. The theory does not take into account collective effects in the ground state. 
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Bose-Einstein condensation (BEC) was discovered (or, 
better, predicted) by Bose ||] and Einstein [|j in 1924- 
25 studying the properties of a gas of non interacting 
bosons in a (infinite) box. BEC has been experimen- 
tally observed only in 1995, in vapors of rubidium ||] 
and sodium Q confined in magnetic traps and cooled. 
Present research on BEC follows similar experimental 
lines, where the potential makes the spectrum discrete 
Q. Theoretical research on BEC in traps is also in 
tumultuous progress , focusing essentially on the role of 
interactions among the atoms. Here we shall instead pay 
attention more to the relation between the discreteness 
and multiplicity of the quantum mechanical spectrum, 
and the features induced by quantum statistics. 

BEC in a box (where the walls do not essentially affect 
the continuity of the spectrum) and BEC in an harmonic 
trap (whose effect is to make the spectrum definitely dis- 
crete), coincide only in some approximate way. First of 
all, in a magnetic trap the number N of confined parti- 
cles as well as the volume V in which they are confined 
are finite, namely we are not in the thermodynamic limit. 
This effect has been discussed in different papers 1^ , Q , 
[ p^ , PI and shown to have limited phenomenological 
relevance, as in the experiments the number of parti- 
cles is 10^ 10^. There is however another reason why 
the thermodynamic limit is not reached - and cannot be 
reached - in magnetic traps, on which we focus our at- 
tention here: the energy spectrum. Quantizing in a box, 
the difference between contiguous energy levels goes to 
zero at least as V~^^^ as volume increases. On the other 
hand, in order to obtain the thermodynamic limit in an 
harmonic potential we should have to require that the 
oscillator frequency - i.e. the spacing of the spectrum 
- goes to zero as N'"^/^. One needs, in other words, 
that the spectrum may be considered continuous. This 
is far from the experimental situation, where typically 
hv/kB ~ IQnK (for v « Hz), to be compared with 
the estimated temperature Tc ~ lO^nK. 

We shall attempt here to propose a better framework 



within which experimental results are consistently de- 
scribed, trying to represent more properly BEC experi- 
mental situation {N finite, energy discrete). Even though 
keeping spectral structure into account may be an impor- 
tant step forward with respect to the continuum theory, 
many body effects are still not considered. The results 
are thus correct up to the detailed structure e.g. of the 
ground state, as it could be described by a collective 
perturbation in the space of free bosons, as in Gross- 
Pitaevsky approach [^,0. 

Ever since 1907 Einstein himself jl^ had pointed out 
how classical results on equipartition of energy are no 
longer valid for T < hv/ks for the quantized harmonic 
oscillator in the frame of Boltzmann statistics. This has 
been so far disregarded, in spite also of the change in 
statistics, as T < hv/ks is below the range of exper- 
iments. Several other authors have on the other hand 
discussed BEC in the context of the present experimen- 
tal situation using equipartition together with the contin- 
uum statistics (eq. (||) below) [|l3|, @, [|l5|, finding that 
the temperature scale of the process is oc N^^^. Within 
the approach presented here this scaling law, essential for 
the coherence of the scheme, is recovered, among other 
novel results. The crucial difference is however that it is 
obtained by re-deriving the entire theory of bosons with 
discrete spectrum, based only on the Gibbs and Bose pos- 
tulates. The central results are that a new equation (eq. 
(^ below) replaces the continuum distribution function, 
equipartition is valid only for T > Tc, and a new critical 
temperature is found. 

Our motivation for this construction is the recent defi- 
nition of bosons as coherent states of the algebra su{l, 1), 
given in |l^]. In that paper a possible generalization of 
the Bose definition of bosons has been considered. Here 
this is used only to define the undegenerate states of the 
harmonic oscillator as the limit for 5 ^ 1 of the general 
formula ( eq. (|l|) below). 

We begin by summarizing the most relevant results 
that a thorough numerical analysis performed within this 



1 



conceptual scheme leads to, leaving for the second part 
of the paper some clarifying details. There turn out to be 
seven essential features when comparing the continuum 
and discrete cases: 

1) The formula (eq. (|^) below) for the occupation num- 
ber n(ei) as a function of the single particle discrete en- 
ergy Ci is the same as in the thermodynamic limit at 
zero energy, moreover for discrete spectrum it does not 
decrease exponentially as e ^ cxo but rather it becomes 
identically zero for e > emax (where Cmax has a well de- 
fined value much less than the total energy E and goes to 
infinity only in the continuum limit). Measurable effects 
are related with the product of n{ei) with the multiplic- 
ity of energy levels m(ei). Fig. 1 displays the product 
n(ei) m(ei) vs. i. 
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FIG. 1. n{ti)m{ti) vs. i for = 10^ bosons in an isotropic har- 
monic trap for different values of T (Tg < T.^ < ) in proximity 
of Tc, showing the appearance of no ^ 0. 

2) The Lagrange multiplier /?', that in the continuum 
is identically equal to \/{kBT) - and oc {E/N)~^ due to 
equipartition - is now for T > Tc /3' oc {E/N)^'^, yet the 
equipartition theorem is still satisfied (constant specific 
heat). In other words, the Lagrange multiplier turns out 
to be a more complicated function of T than oc T~^. 

3) Temperature, defined - without resorting to 
equipartition - as the derivative of energy with respect to 
entropy, allows us to describe entirely the thermodynam- 
ics of the system, assumed in global thermal equilibrium 
(as done in |]l3|) because of the efficiency of evaporative 
cooling. For T > Tc, effects related to the spectrum 
discreteness appear to be irrelevant, as they are in deter- 
mining Tc as a function of N , . Another critical 
temperature Td is however found (see Fig. 2) defined by 
the ground state filling uq/N crossing the value 1/2. 

As Fig. 2 shows, C is constant, equal to Sfc^ (as it 
should, since the system does not condense and behaves 
as a classical gas) for T > T^, it is almost zero (as en- 
tropy and energy decrease with quite different slopes be- 
cause only a few atoms progressively migrate from the 
warmer tail of the spectrum into the fundamental level) 
for Td < T < Tc, and exhibits a peak for T = Td, where 



condensation becomes a global effect. Note that, after 
the spike at Td, the specific heat goes to zero as T goes 
to zero, as required by Nernst's theorem. 
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FIG. 2. Specific heat C vs. T in units kB = hu = 1 ior N = 10^ 
bosons in an isotropic harmonic trap, as predicted by discrete quan- 
tum statistics. 

4) The chemical potential is no longer equal to the La- 
grange multiplier a' and does not appear to present ma- 
jor differences from the results reported in the literature 
§■ 

5) Repeated numerical analysis for various values of N 
indicates that Tc oc iV^/'^, in agreement with ||T^, [ p5[ , 
whereas one finds Td oc N^/^/\nN. 

6) As explicitly shown in Fig. 1, the high energy tail 
of the velocity distribution is not of Maxwell-Boltzmann 
type and it cannot be straightforwardly used to deter- 
mine the experiment temperature |13|. 
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FIG. 3. The fraction of condensate no/N vs. Af'^/^ {E/Nhu-3/2) 
for N = 10* (• ■ ■), 10"^ (- - -), 10« (— ). 

7) Thermo-dynamical temperature, despite its concep- 
tual relevance, does not appear to be the physical param- 
eter that best describes the process of condensation: it 
is difficult to measure in the typical experimental set-up, 
and as one does not have equipartition its intuitive mean- 
ing of "average energy per particle " is lost. We prefer 
then to describe the onset of condensation by exhibit- 
ing the behavior of no/N vs. the true (dimensionless) 
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"energy per particle" E/Nhv (Fig. 3, where in abscissa 
we actually reported N~^/^ (^E/Nhv — 3/2) to point out 
the dependence on N^^^), much easier to compare with 
experiments. Our results are consistent with the experi- 
mental findings in ref. Fig. 4, of which Fig. 2 shows 
the derivative, provides E/N as a function of T . This 
curve can be inverted giving T = T(E/N), a relation 
which replaces equipartition (where it would be linear) 
and allows us to evaluate temperature dependent quan- 
tities as functions of the experimentally more accessible 
internal energy. 
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FIG. 4. Specific internal energy E/N vs. T, for = 10® (units as 
in Fig. 2). 

We shall now discuss in some mathematical detail how 
the above features come into play. The traditional text- 
book way (see e.g. [0) of deriving the boson gas equilib- 
rium distribution is the following. Starting from Gibbs' 
assumption, in order to derive the combinatorics many 
single particle states (levels) of energy practically equal 
are collected together in cells (an assumption justified a 
posteriori by the thermodynamic limit) and the parti- 
cle distribution is straightforwardly found to be, up to 
normalization, given by 



w{n,} = n 



r(nj -t- g) 

r(n, + i)r(5) 



where g is the multiplicity, that with no loss of gener- 
ality we assume equal for all cells. The equilibrium dis- 
tribution, identified with that of maximum probability, 
is then obtained by finding the set {rii} that maximizes 
(the logarithm of) VF{ni}, constrained by the two con- 
ditions Ui = N and eiUi = E. One obtains the 
exact equation 



i'i'rii + g) - ipirii + 1) = a + /Sei 



(1) 



where a = a{g), dimensionless, and (3 = P{g), with the 
physical dimension of an inverse energy, are the two La- 
grange multipliers, clearly depending on the value of g, 
introduced in order to implement the constraints, while 
t/j is the digamma function pi]. 



For large argument ■0 coincides with the logarithm, 
thus for Ui ^ 1 and g 3> 1 it is easy to obtain (identify- 
ing n{ei) with the finite density rii/g) 



n{ei 



1 



(2) 



The latter is the Bose-Einstein distribution we are accus- 
tomed to, so well established that it is sometimes consid- 
ered as an alternative definition of bosons. The essential 
assumption of this derivation is of course that the energy 
spectrum is continuous, hypothesis that does not hold in 
the presence of an harmonic potential. 

Eq. (|l|) has more general validity in that it holds true 
for any value of g and hence also for the discrete case. For 
example, in the case of photons in a laser beam, where 
degeneracy of states is two because of helicity, it gives 

@ 

ni = 1 for hi < 1 and n-i = for hi > 1 . (3) 

iii 

For a physical (3-dimensional) quantum harmonic os- 
cillator, g should be set equal to 1. Indeed, even in the 
"isotropic" case, that we shall consider as an example, 
the experimental frequencies are such that Vi — Vj t'i, 
but Vi is never exactly equal to Vj] we have then a com- 
puting degeneracy only, because statistics remains un- 
degenerate. However, setting <? = 1 leads for ([^) to a 
perfectly acceptable but void equation (left hand side 
identically zero, implying q;(1) = (3{1) = 0). On the 
other hand: i) eq. (|l|) is defined for g in the whole com- 
plex plane except for negative integers, and ii) in ref. 
an alternative definition of bosons has been proposed in 
terms of coherent states of the discrete series of SU (1,1), 
where the highest weight of the representation is g/2, 
with g any strictly positive real number. In this latter 
context g is not restricted to the integer values imposed 
by the combinatorial interpretation and we can define by 
analytical continuation for g — > 1 the equation for unde- 
generate levels that we need, rii = n{ei) results thus to 
be the solution to the equation 



i}'{ni + 1) = a' + P'e, = h'. 



(4) 



where a' = lim a(g)/{g — 1) and /?' = lim (3{g)/{g — 1) 

s— 1 3-»i 
are the new Lagrange multipliers to be determined by 

the constraints of fixed N and iJ, and ip' is the trigamma 

function [Q, derivative of ip- Eq.(^ is one-to-one and 

can be inverted, giving us the correct formula for discrete 

spectrum, to be used instead of eq. (0) : 



and 



= [,p']~^{h'i) -1 for h'^<n'^/6 



= for h'^ > , 



(5) 



[ip'] ^ denoting the inverse function of ■0'- tt'^ / 6 = ip' (1) . 
For /i' — » eq. (||) has the same asymptotic behavior as 
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: Tii ^ h'~^ — ^. The two formulas (^) and (||) ex- 
hibit a similar structure, characteristic of all values of g, 
with a pole (with residue 1) at argument zero and a zero 
for large values of the argument (compare eq. (^). The 
discrete description, however, does not require an ad hoc 
delta function at the ground state, which is instead dis- 
cussed on the same footing as all the others. In fact one 
may assign first the values of N and E; from these the 
two Lagrange multipliers a' — a'{N,E), (3' — I3'{N,E) 
are then obtained from the constraints using for Ui the 
defining equation (|^); then, inserting them in (^ itself, 
one gets rn = ni{a' , (3') = ni{N,E) (see Figs. 1 and 3). 

As condensation derives from the existence of the pole 
combined with the feature that multiplicity increases 
quadratically for the excited states, the continuum for- 
mula (^ is correct in predicting the collapse of the atoms 
in the fundamental level, but it is unable to describe 
what happens for T < Tc- The subtlety here is that 
the very concept of temperature, as defined in the the- 
ory of gases, cannot be straightforwardly extended to the 
discrete spectrum case. 

Numerical calculations performed for a physical 3- 
dimensional "isotropic " quantum harmonic oscillator 
with N ranging from 10^ to 10^^, made for simplicity 
approximating [ip'] 



for hi < 7r^/6 and n 

h, TT^ 



^ by fixing only pole and zero, 

for h^ > TrVe , 



show indeed that, for E/N large, 7r^/6 



a oc 



while /?' cx {E/N) . The internal energy is therefore 
not proportional to the inverse Lagrange multiplier /3' 
but to the more complex expression (tt^/G — a')/3' ^. 
Temperature can then no longer be obtained, as in the 
continuum, by T = {ksfi') and equipartition of en- 
ergy must be demonstrated. The only way to proceed 
is to adopt the basic definition of temperature based on 
entropy. In order to do so, we use Shannon's definition 
of entropy in terms of probabilities 



{E/N) 



-3 



where we set pi = Ui/N, recalling that the ni{N, i?)'s are 
now known. At fixed N we thus straightforwardly obtain 
S — S{E). Simple numerical operations finally give the 
temperature 

T = T{E) = 



dE 



N 



Once more the function is one-to-one, and from T ~ 
T{E) we can obtain E = E{T) (see Fig. 4), whence the 

is obtained by (numerical) 



specific heat C = 



1 

N 



dT 



N 



differentiation (see Fig. 2). Also fi, which is no longer 
related in a simple way to a', can be calculated in anal- 
^dE] 

and all the results presented 

s 



ogous way, fi 



dN 

above are thus obtained 



It should be pointed out that the theory presented 
leads not only to a set of specific behaviors, but it de- 
termines as well the measurable quantities: specific heat, 
filling order parameter, chemical potential, etc., all of 
which can be experimentally measured and are related 
to the characteristics of the experiment {N, E, v). To 
this latter effect, it should be pointed out that in the ex- 
periments the harmonic potential is usually anisotropic: 
extension of our results to the corresponding physical pa- 
rameters is straightforward. Whether or not the features 
predicted theoretically are correct or not can therefore be 
tested by direct measurements of sufficiently high sensi- 
tivity. This feedback from the laboratory is crucial before 
further refinements of our analysis, such as e.g. inclusion 
of many body effects, is carried over. 
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